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I. INTRODUCTION 

It is very difficult to describe full electron correlations 
due to a number of complex phenomena related to the 
quantum character of the electron. It is hence appro- 
priate to approach the full description of electron cor- 
relations iteratively in several stages. The first one is 
the static mean-field approximation of the Hartree type. 
Such a mean-field theory completely neglects both charge 
and spin fiuctuations and reduces the interacting system 
to a Fermi gas with renormalized, self-consistently de- 
termined, densities. Such a simplification may deliver 
reliable results only for macroscopic static quantities in 
] the weak-coupling limit. Electron correlations in this ap- 
proximation have no impact on dynamical and transport 
properties. 

The next level in a comprehensive modeling of strong 
electron correlations are models and approximations al- 
lowing for charge fluctuations. In these models the spin 
of electrons does not play a significant role and electrons 
' are subjected only to a potential scattering. It means 
that electrons are only scattered on fluctuations of the 
atomic potential in the lattice. The potential fiuctua- 
tions are caused by impurities that may be distributed in 
the crystal either regularly or randomly. The paradigm 
for the former situation is the Falicov-Kimball model^ 
(FKM) and for the latter the Anderson model of dis- 
ordered electrons (DAM).^ Unlike the static mean-field 
approximations, the models with a potential scattering 
lead to quantum dynamical effects and are applicable to 
the entire range of the interaction strength (variance of 
the potential fiuctuations) . The potential scattering does 
consequently affect spectral and transport properties of 
the system. 

A common feature of the models with potential scat- 
terings only is that energy is conserved during scattering 
events and need not be treated as a dynamical variable. 
Each energy, however, is renormalized in a different man- 
ner and hence the energy (frequency) is used as an ex- 
ternal label. Conservation of energy in scattering events 
is a significant simplification in the description of elec- 
tron correlations. It allows for an exact solution in the 



limit of infinite spatial dimensions (dynamical mean-field 
theory), where the effect of strong potential fluctuations 
may be studied without uncontrolled approximations!^^— 

The two models, FKM and DAM, are standardly used 
for different purposes. The former one is aimed at a de- 
scription of quantum fluctuations caused by electron cor- 
relations in thermally equilibrated states, while the latter 
one was introduced so that a response of a disordered 
electron gas to weak electromagnetic non-equilibrium 
perturbations can be estimated in a controlled way. Both 
the models have served well their original purposes. 
The Falicov-Kimball model has been successfully applied 
to a simplified description of correlation-induced metal- 
insulator^'^ and valence-change^ transitions in rare-earth 
compounds, or atoms in optical lattices. The dis- 
ordered Anderson model has been used to describe the 
spectral and transport properties of metallic alloys^- and 
vanishing of diffusion, called Anderson localization^^ 
There have been efforts to describe the combined ef- 
fect of electron correlations and randomness in the dis- 
ordered Falicov-Kimball model^ or Anderson localiza- 
tion in FKMJ^ Only a few attempts have, however, 
been made in the calculation of the response of FKM 
to non-equilibrium perturbations beyond the mean-field 
approach. In particular, it is little known about the 
electrical conductivity of FKM beyond the mean-field, 
Drude contribution^^ 

It is the aim of this paper to fill up this gap and to 
propose a systematic way how to calculate vertex correc- 
tions to the Drude (mean-field) electrical conductivity in 
models with elastic scatterings only, that is FKM, DAM 
or a disordered FKM. The method we use is an expansion 
around the mean-field solution obtained from the asymp- 
totic limit to high spatial dimensions. We calculate the 
leading-order vertex correction in high spatial dimensions 
being of order while the Drude conductivity is of 

order 1/d. We demonstrate that the vertex corrections 
have a universal behavior for all models of elastically scat- 
tered electrons and are always negative, independent of 
whether they are caused by an external random poten- 
tial (quenched randomness) or by static, thermally equi- 
librated electron correlations (annealed randomness). 
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II. ELECTRICAL CONDUCTIVITY OF 
ELASTICALLY SCATTERED ELECTRONS 

Elastic scatterings of electrons are caused by either 
internal or external fluctuations of atomic potentials of 
frozen ions forming a crystalline lattice. That is, interac- 
tions of electrons of the same sort are excluded in models 
with elastic scatterings. They are actually forbidden in 
spinless models with locally interacting fermions. These 
models hence contain either more than one type of elec- 
trons or an externally governed distribution of atomic 
potentials. We consider only the elementary version of 
such models and take into account only two types of elec- 
trons, extended and localized ones. We further assume 
homogeneity in the distribution of the localized electrons 
and hence the non-interacting part of the Hamiltonian 
describing such a situation reads 



UJ 



q, UJ 



ifo-E^(k)ct(k)c(k) + £;^5]///, 



(la) 



We are interested in dynamical properties of the delo- 
calized electrons induced by fluctuations of the atomic 
potential the extended electrons feel. To this purpose we 
introduce an interacting term 



(lb) 



where q — N~'^ c(k) exp{— «k • R^}. We denoted Vi 
atomic levels of the ion situated in the elementary cell 
centered around the lattice vector R^, and U is the inter- 
action strength between the extended and localized elec- 
trons. We generally assume that the atomic potential 
Vi is a random variable with a static site-independent 
probability distribution of its values determined exter- 
nally. If t/ = the full Hamihonian H = Hq + Hj 
is that of the Anderson model with disordered electrons 
and if Vi = 0, the full Hamiltonian describes the Falicov- 
Kimball model. We hence see that the most general ele- 
mentary Hamiltonian for elastically scattered electrons is 
just a disordered FKM. Actually, both the contributions 
to the interacting Hamiltonian from Eq. (jlbp introduce 
randomness into the distribution of atomic potentials the 
extended electrons feel. Potential Vi represents a static 
(quenched) randomness and the interaction U a dynam- 
ical (annealed) one. Both contributions can be treated 
on the same footing. 

We are not interested in this paper in equilibrium 
thermodynamic properties of FKM but rather in its dy- 
namical behavior and particularly in the static, opti- 
cal conductivity. The electrical conductivity, although 
static, is nevertheless a dynamical property, since we need 
at least two different energies (small imaginary parts) to 
determine it. The Kubo formula for the diagonal part of 
the electrical conductivity in models with only elastically 
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FIG. 1. Translationally invariant (averaged) two-particle 
Green function in momentum space in the notation used in 
the Kubo formula for electrical conductivity. 



scattered electrons can be written asii 
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where Va(k) — de{]i)/dka is the group velocity in the 
direction a and the averaged (translationally invariant) 
two-particle Green functions are defined 



Gi^i^{uj,uj';q) = G\;^,{Lo + tO 
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,w' + iO+;q), 



We denoted oj = E — the energy measured from the 



Fermi level and f{E) = 1/(1 



) is the Fermi 



function. See Fig. [l]for the way the variables in the two- 
particle Green function are used in this paper. 

The expression for the electrical conductivity simpli- 
fies if we resort to zero temperature. Then the integrals 
over frequencies can be performed explicitly and we ob- 
tain a simple formula 

v^{\.W{\^)[Gt5-'^G^^] (3) 



kk' 

with the values of the two-particle Green functions at 
the Fermi energy. We used an abbreviation G^^, = 
Gi^«(0,0;O). 

The two-particle Green function G*-^^ carries informa- 
tion on both the uncorrelated and correlated motion of 
two electrons. Only the latter one is the genuine two- 
particle quantity. It is our task to identify this contribu- 
tion to the electrical conductivity. To do so, we introduce 
the two-particle vertex F defined from an equation 



'-'kk' — Gk Gk [^(k — k') + Fi^i^ G^/Gk'] 



(4) 



where we again used a notation G^ = G^(0). With the 
aid of the vertex function we can decompose the conduc- 
tivity tensor into two parts 



.(0) 



(5) 
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where 



.(0) 
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TTiV ^ 
k 



«„(k)n3G«(k: 
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is the standard one-electron or Drude conductivity at 
zero temperature. The genuine two-particle contribution 
is called a vertex correction and is proportional to the ap- 
propriate matrix element of the two-particle vertex that 
at zero temperature reads 

Aa„„ = ^ E^"««"(k') {|G^r Ari^« \G^,\' 
[G^y AT^^ (G^,)']} . (7) 
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It is not the full two-particle vertex F that is important 
for the electrical conductivity, but only its odd part AF. 
That is, only the part of the vertex function being on 
bipartite lattices an odd function in fermionic momenta 
k and k', contributes to the electrical conductivity. 



III. EQUILIBRIUM MEAN-FIELD 
THERMODYNAMICS 

The Kubo formula for the electrical conductivity, 
Eq. was derived within the linear- response theory 
and hence the Green functions entering this formula are 
the equilibrium ones. To estimate quantitatively the ver- 
tex corrections to the electrical conductivity we need 
to know the equilibrium thermodynamics of FKM. To 
reduce the impact of uncontrolled approximations we 
should at best know the exact equilibrium grand poten- 
tial. An exact solution to FKM is known in the limit of 
infinite spatial dimensions (mean-field limit) with which 
we start up. The equilibrium thermodynamics of the 
disordered FKM in d = oo was analyzed in Ref. The 
functional of the averaged grand potential was found to 
be represented via a set of complex variational parame- 
ters Gn and E„, where the index n corresponds to the 
nth fermionic Matsubara frequency 



+ (ln[l + G„(E„-y)]), 

-{\n[l + e^p{/3{fi-Ej-£v}]),., ■ (8) 

Symbol ( stands for averaging over the distribution 
of the random potential Vi. The shift of the /-electron 
atomic level Sy is determined via the same complex num- 
bers Gn and E„ 
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and depends on the configuration of the random atomic 
potential V . 

The equilibrium thermodynamics is obtained as a sta- 
tionarity point with respect to small variations of com- 
plex numbers S„ and G„ of the averaged grand potential 
(fj)^^ from Eq. ([5]). Vanishing of variations of the former 
and the latter parameters lead to a couple of equations 
for each Matsubara frequency za;„ 



Gn 



dep(e) 



l~fiEf+£v) , f{Ef+Sv) 



(10a) 



l + G„(S„-y) l + G„(E„-F-C/) 
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The first equation states that in equilibrium G„ is the 
local element of the one-electron thermal Green function 
with a self-energy E„. The second equation determines 
the value of the equilibrium self-energy. These equations 
of thermal equilibrium must be completed with an equa- 
tion determining the chemical potential p from the total 
electron density n. This equation then is 

oo 



{fiEf+£v))a. 

^/(c.)5G«(^-E^(..)) 



(11) 



where G«(z) = iV-i G^(k, z). Equations 

fully determine the equilibrium thermodynamics for a 

given temperature T and a total particle density n. 

Only one-particle equilibrium functions can be di- 
rectly calculated from the grand potential {^)av- To de- 
rive higher-order correlation functions we have to slightly 
perturb equilibrium and look at the corresponding re- 
sponse functions. For the electrical conductivity we need 
to know two-particle vertex F. The only consistent way 
to derive a two-particle vertex within the mean-field the- 
ory is to keep the non-equilibrium perturbation local. — 
The resulting vertex remains local and is only frequency 
(energy) dependent. 

The equilibrium two-particle vertex has generally 
three independent (Matsubara) frequencies. In the case 
of FKM investigated here, the resulting two-particle ver- 
tex can have maximally two independent frequencies. 
The two frequencies can, however, be placed in two dif- 
ferent ways. The full local two-particle vertex for FKM 
can be represented as 



F 



MF 



(12) 



The two contributions to the full vertex correspond to 
two ways the electron lines go through the vertex. The 
electron line entering into the vertex at the upper left 
corner goes out either via the upper right corner (vertex 
7) or via the lower left corner (vertex ip). See Fig. [5] for 
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FIG. 2. Graphical representation of Eq. (|12p . Tiie dasiied 
lines within the boxes indicate charge propagation from the 
incoming to the corresponding outgoing line. 



a graphical representation and notation used in Eq. ()12|) . 
Due to energy conservation in this model the incoming 
and outgoing frequencies must equal. Notice that the 
two vertices 7 and never mix up and are completely 
detached in the solution. The former vertex is relevant 
for the transport while the latter one for the thermody- 
namics. It is hence sufficient to take into account only 
vertex 7 for the calculation of the electrical conductivity. 

We can represent vertex 7 via a local Bethe-Salpeter 
equation with a local irreducible vertex h.m,n- We have 



1 ^m^nGmGn 



(13) 



Charge conservation leads to a generalized Ward iden- 
tity that matches a non-equilibrium variation of the self- 
energy <5Em,n with the equilibrium irreducible vertex 



A. 



A, 



1 



where we abbreviated 



[1 _ \-l 



(14) 



1 - f{Ef + fv 



[1 + Gmi^m - V)] [1 + G„(E„ - V))] 

fjEf+Sy) 

[1 + G„(I]™ -V-U)][l + G„(S„ ^V- C/)] 

(15) 

The latter representation was derived with the aid of sta- 
tionarity equations, Eqs. (fTO)) . 

We have derived all the necessary equilibrium quanti- 
ties needed for the calculation of the electrical conductiv- 
ity. We must, however, go beyond the local, mean-field 
limit in order to calculate transport properties. That is, 
wc have to perturb the equilibrium with time-dependent 
non-local excitations. To do it in a controlled manner 
wc use a perturbation expansion around a mean-field. 
This can be consistently achieved within an asymptotic 
expansion in high spatial dimensions. 



IV. EXPANSION AROUND MEAN FIELD - 
VERTEX CORRECTIONS TO THE ELECTRICAL 
CONDUCTIVITY 

An asymptotic expansion around a mean-field solu- 
tion in d = 00 for non-interacting disordered electrons 



was recently derived by one of us in Ref. [191 The only 
ingredients of such an expansion are the local mean-field 
vertex 7 and the the non-local one-electron propagator. 
The relevance and applicability of this expansion goes 
beyond the model of disordered electrons. It can be ap- 
plied to any model where energy is conserved in scat- 
tering events (elastic scatterings) and where only spatial 
fiuctuations matter. The disordered FKM studied here 
falls into this category^ and we can hence use the expan- 
sion concept of Ref. ,19 for it. It is a considerable ad- 
vantage to expand around a mean-field solution, since all 
the effects of scatterings on impurities are included al- 
ready in the local vertex 7. The variance of the potential 
fiuctuations is not a small parameter and the expansion 
terms do not depend on whether the scatterings are due 
to thermally equilibrated or randomly distributed static 
impurities, that is, whether the disorder is annealed or 
quenched. 

We construct the expansion around a mean field as an 
asymptotic series on a hyper-cubic lattice in high spatial 
dimensions. The expansion parameter is the off-diagonal 
one-electron propagator from the mean-field theory. We 
define 



G(k,C) 



1 



C-£(k) 



depje) 



(16) 



where we denoted = z — S(z) and the local self-energy 
S(z) is that of the mean-field solution. The off-diagonal 
two particle bubble is a convolution of the off-diagonal 
one-electron propagators. We hence define 



x(C,C';q) = ^E^(i^'C)G(k-Hq,c' 



(17) 



The frequency indices are external parameters and we 
suppress them when they are not necessary to specify a 
particular type of the one or two-electron propagators. 

The asymptotic limit of the full two-particle vertex 
in high spatial dimensions contains beyond the local 
mean-field vertex 7 also non-local contributions from the 
electron- hole and electron-electron ladders. It can be 
represented as follows 



kk'lqj 



7 



1+7 



x(q) 



x(Q) 



1 - 7x(q) 1 " 7x(Q) 



(18) 



where we denoted Q = q + k + k' the momentum con- 
served in the electron-electron channel. Notice that the 
contribution from the electron- hole channel with x(q) is 
part of the two-particle vertex from the Coherent Po- 
tential Approximation (CPA) and can be derived from 
a Velicky-Ward identity.^" The two-particle vertex from 
CPA does not carry the full 1/d correction to the local 
vertex and moreover it is not electron-hole symmetric 
on the two-particle levelJ^ A consistent extension of the 
local mean-field two-particle vertex must contain both 
non-local contributions from the electron-hole and the 
electron-electron channels as given in Eq. ([T5)) . 
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The contribution from multiple scatterings in the 
electron-electron channel to the asymptotic two-particle 
vertex is important in particular in the calculation of the 
vertex corrections to the mean-field electrical conductiv- 
ity. The CPA vertex, as notoriously known, does not 
generate vertex corrections, unless we introduce odd dis- 
persion relations in multi-orbital models. It is only the 
second term in the parentheses on the right-hand side of 
Eq. (flSl) . that contributes to the averaged conductivity 
and the we can identify 



Arkk'(q)=7' 



X(k + + q) 
1 - 7x(k k' -f q) 



(19) 



from the formula for the vertex corrections to the electri- 
cal conductivity, Eq. ([T]). Inserting representation from 
Eq. (HH) into Eq. ^ we obtain 



kk' 
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2 1^,2 V^m^'^i^ + ^'Wi^') 



G^'f I 



RR\^ (r<R\^ (k + k')i;a(k') , ^^.2 



1 — ^RR-y^RR{\i + k' 
It is a general formula for the leading-order corrections 



to the mean-field (Drude) conductivity and can be ap- 
plied in any dimension. Vertex AF contains the so-called 
Cooper pole being the image of the diffusion pole con- 
tained in the CPA vertex after an appropriate electron- 
hole transformation. The Cooper pole is responsible for 
the so-called weak-localization corrections in the Ander- 
son model of disordered electrons.— These corrections are 
here identical with those determined by Eq. (pUj) . They 
are negative and diverge in spatial dimensions d < 2. The 
mean-field corrections to the electrical conductivity from 
Eq. (|20|1 can effectively be applied only in dimensions 
d > 3 and not too close to the band edges. 

We can further simplify Eq. (|20p in that we get rid 
of the local mean-field vertex 7. We utilize the Ward 
identity from Eq. (fT4|) connecting the one-electron self- 
energy S and the two-particle irreducible vertex A. Using 
the stationarity equation, Eq. (jlOal) . we easily obtain 



A. 



A, 



1 



G'+ x++(0) 



X+ 



(0) 



(21a) 
(21b) 



where we denoted = dT,+ {uj)/duj\u:=o, G'_^ — 

dG+{uj)/duj\^=o, and Z+ = ^^/(S; - 1). With these 
(20) relations and the local Betlic-Salpeter equation, Eq. ([T3|. 
we can rewrite Eq. (1^ to 



A MF _ ^ 
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\a(k)x+- 
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k')«a(k')|G^,|' 


(x+-(o) 


-\{G+)f^ 


)(x+ 


_(0)-x+-(k + k')) 
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Zl (Gg)%^(k)x++(k + {G^.f 
(x++(0) - Z+(G+)2) (x++(0) - Z+x++(k + k')) 



(22) 



where we used an abbreviated notation 



{Gi 



dep{e) 



_^ [EF~e~J:±± iO+] 
dtp{e) 



-00 [Ef -e-^±± iO+] 



2 ■ 



(23a) 
(23b) 



Representation of the vertex corrections to the 

electrical conductivity does not explicitly contain the 
strength of elastic scatterings in the model. This strength 
is beyond the self-energy of the one-electron propaga- 
tors implicitly comprised in the spatial fiuctuations of 
the two-particle bubble x(q)- A singular structure of the 
integrand in momentum representation of the vertex cor- 
rections to the Drude conductivity becomes transparent 
in Eq. (HH). 



V. VERTEX CORRECTIONS TO THE 
ELECTRICAL CONDUCTIVITY FROM HIGH 
SPATIAL DIMENSIONS 

To evaluate the vertex corrections to the mean-field 
electrical conductivity, Eq. ([20]), we resort to high spa- 
tial dimensions where we can explicitly perform the con- 
volutions over momenta. We are interested only in the 
leading-order contributions in the inverse spatial dimen- 
sion 1/d. The conductivity vanishes in the mean-field 
limit, d = (X), where only local quantities survive. The 
actual mean-field conductivity is due to the velocity in a 
particular direction proportional to 1 /d and hence can be 
treated only asymptotically for d — > 00. On a hyper-cubic 
lattice we have VaCk) = de(k)/dka = td~^^^ ainka- In- 
serting this result in the mean- field conductivity, Eq. (jG]) , 
we obtain the Drude conductivity on a hyper-cubic d- 
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dimensional lattice to be 



(24) 



To evaluate the vertex corrections we represent the 
one and two-particle propagators so that we can sep- 
arate the Cartesian components of momenta. We use 
the following integral representation for the one-electron 
propagator 

G(k, C±) = -i du e±™^± TT exp{±^ cos K} , 
Jo V d 

(25) 

where we assumed that 5C+ > and 5^- < in or- 
der to keep the integrals convergent. Here ki, is the I'th 
Cartesian component of momentum k on a d-dimensional 
hyper-cubic lattice. 

The two-particle bubble can be represented in a sim- 
ilar way. Performing the integration over momenta we 
obtain in the leading asymptotic order for d — > oo 



x(q;C,C') = - / du I dv e-'-e 

Jo Jo 



exp{ 



uvt'^ 



cosg^} . (26) 



For simplicity we assumed that both complex energies C 
and ^' have positive imaginary parts. A generalization 
to different imaginary parts is straightforward. 

The above integral representations suffice to evaluate 
the leading order of the vertex corrections. We realize 
that the denominator in the representation of the vertex 
corrections, Eq. (I^Ul) . does not contribute in the leading 
order for d oo. We denote 

= (G^O'-oWx^'^Ck + kO^.lk') (G^,)' . 



kk' 
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We use the integral representation separating the Carte- 
sian components of momenta. We then have 
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J«« = - / da db da' db' / du 
d Jo Jq Jo Jo Jo 



exp 



X sin ka sin k'^ JJ^ exp < 



i{a + b)t 
Vd 



cos fcjy 



+ 



i{a' + b')t 
Vd 



cos k' 



uvt 



cos(fcy + k'^) 



The leading non- vanishing contribution from the summa- 
tion over momenta is proportional to 1 / d. Performing the 



calculation we obtain 



Jfi* = ^ / da \ db I da' \ db' I du 
°d Jq Jo Jq Jo Jq 



dv e'('^+'')C+ 



exp 



(a + b)H^ ^ ^^ia'+b'K+ 



X exp 



X e™^+ exp 



vH-" 



(28) 



Using representation 
form 



J, 



RR 
aa 



we can rewrite J^-^ in a simple 



{G+G+f{Gl){Gl) (29a) 



where each Green function G+ stands for the selection of 
the imaginary part of the propagator, the retarded one. 
In the case of the retarded and advanced propagators we 
then have 



J. 



RA 



{G+G^f{Gl){Gl) 



(29b) 



We use these results to represent the vertex correction 
to the electrical conductivity in high spatial dimensions. 
We obtain 



2+4 



ArfCr 
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(G^)Si+]} . (30) 



With the aid of the Ward identity, Eq. (PT|) , we can fur- 
ther simplify the expression for the leading contribution 
to the electrical conductivity beyond the Drude formula 
in d spatial dimensions 
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{Gl)-Z+{G+f 



(31) 



The sign of the vertex correction is negative. Al- 
though it is only the leading asymptotic term, it de- 
termines the overall sign of the vertex correction. The 
neglected higher-order terms must be summed to infinite 
order so that the denominators on the right-hand side of 
Eq. (f22|) are recovered. The higher-order terms then do 
not change the sign of the leading-order vertex correction 
to the electrical conductivity. 



VI. RESULTS 



To reach numerical values for the vertex corrections 
to the electrical conductivity we need the stationarity 
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equation for the self-energy, Eq. (jlObp . This equation 
significantly simplifies in the pure model at half filling 
with a symmetry between c and / electrons, it means 
if He = 1/2, Uf = 1/2, Ef = 0, and Ep = U/2. We 
then obtain 3?E(0) = U/2 and the solution for the local 
functions can be expressed in terms of a single parameter 



depje) 



> 



It is easy to find from Eq. (|10bl) that 



3S+ = - 



3G+ 



7+- 




1++ 



(32a) 
(32b) 
(32c) 



The stationarity equation for the parameter x reads 
dep(e) _ / 1 \ 



1 = 



X 



4 



(33) 



where we again abbreviated the integral over energy 
weighted by the density of energy states by angular 
brackets. 

The two-electron bubble for zero momentum can be 
represented as follows 



{G\) = . 



1 - — X + 2 ( 1 - — X I X2 



(34a) 



where we introduced variance of non-local fiuctuations 



X2 



dep{e) 



The Drude conductivity then is 



and the vertex correction reads 
Ada 



(34b) 



(35) 



e 

X <! U^x - 



C/4 

4a;2 



C/2 



1 —X + 21 —X X2 



C/2 



1 - —X I X2 



(36) 



We set the energy scale t = 1. Note that < x < A/U^ 
for [/2 < 4(e2). The upper limit on the interaction 
strength U in the above equations is imposed by the 
metal-insulator transition at which x — and the density 
of the extended electrons at the Fermi level vanishes. 



Pf 



U 



FIG. 3. (Color online) Drude conductivity vs. density of 
states for the electron-hole symmetric Falicov-Kimball model 
with semi-elliptic density of states and d = 3. 



The Drude conductivity in the weak-coupling limit 
{U — > 0) diverges, that is the resistivity vanishes as 
it should be for the Fermi gas without impurity scat- 
terings. The two parameters x and X2 behave in the 
weak-coupling limit x ^ A/U"^ — tt"^ p{QY and X2 ^ 
2/U^TT^p{0 f - 1 and hence 



CTO 



^ IP 



(37) 



The vertex correction from Eq. p6p remains finite in the 
weak-coupling limit. To derive an explicit expression for 
it we had to expand the asymptotic solution for the pa- 
rameters X and X2 up to the third order in d^^. 

For the explicit calculation we used the semi-elliptic 
density of states p{uj) = 2/tt\/1 — uj^. We choose the di- 
mensionality parameter d = 3 in our calculations. The 
Drude conductivity is plotted in Fig. [31 It is compared to 
the density of states (DOS) of the mobile electrons at the 
Fermi energy. Coulomb interaction decreases both the 
DOS and the conductivity down to the metal-insulator 
transition, where the self-energy diverges. The vertex 
correction Act from Eq. (p6)) is plotted in Fig. |4l The 
modulus of the vertex correction is not a monotonic func- 
tion and it reaches maximum at about [/,„ ~ 0.82. The 
vertex correction is negative but is much smaller than 
the Drude conductivity. Their ratio is plotted in Fig. m 
Only close to the metal-insulator transition the vertex 
correction is of order of the mean-field conductivity. 

The Falicov-Kimball model away from half filling is 
more complicated. First, we need two parameters to de- 
scribe the local quantities. They are 



= (|G^ 



dep{e) 



{Ep - 3?S - ey + 
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-0.018 I 1 ' 1 1 
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U 

FIG. 4. Vertex correction to the conductivity for the electron- 
hole symmetric case. 



0.7 




FIG. 5. Ratio r — \Aa\/ao for the electron-hole symmetric 
case. 

and 




dep(e)e 



where we subtract the Fermi energy and the self-energy 
from their values in the electron-hole symmetric case 
rif ~ 1/2, Tic = 1/2. Explicit formulas are more involved 
and we do not present them here. Second, as discussed 
in Ref. 22, a nontrivial solution at zero temperature ex- 
ists only for 1/2 < n < 3/2. Outside this region the 
/-electron energy level is either empty (n < 1/2) or fully 
filled (n > 3/2). We plotted the Drude conductivity and 
DOS in Fig. [5] for a total filling n = 0.7. The density 
of states at the Fermi energy is almost constant within 
the interval < U < 0.82, within which the extended 
electrons scatter on /-electrons. The extended electrons 
go over to a Fermi gas in both limiting values of inter- 
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FIG. 6. (Color online) Drude conductivity vs. density of 
states for filling n = 0.7. 
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-0.01 I ' 1 1 ' 1 ' 1 J 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

U 

FIG. 7. Vertex correction to the conductivity for filling n = 
0.7. 

action, hence the Drude conductivity diverges at both 
ends. The asymptotics is, however, different at the two 
ends. The asymmetry between the two limiting interac- 
tion strengths is demonstrated in Fig. [71 where we plot- 
ted the vertex correction to the Drude conductivity. It 
is two orders smaller than the mean-field one. This is 
clearly seen from their ratio plotted in Fig. [S] The ratio 
is no longer a monotonically increasing function as in the 
electron-hole symmetric case, but reaches maximum at 
Um « 0.48. 



VII. CONCLUSIONS 

We studied effects of elastic scatterings of mobile elec- 
trons on either thermally equilibrated or frozen, ran- 
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0.006 




FIG. 8. Ratio r = |Act|/o-o for filling n = 0.7. 



domly distributed static impurities. We concentrated 
on quantum coherence effects due to correlated back- 
scatterings and its impact on the electrical conductiv- 
ity. We calculated vertex corrections to the mean-field 
(Drude) conductivity. We used a systematic expansion 
around the mean-field solution via the asymptotic limit 
to high spatial dimensions. Our principal finding is that 
elastic scatterings always lead to diminution of the Drude 
conductivity. That is, vertex corrections due to elastic 
scatterings have negative sign. The sign of the vertex 
correction is determined by its leading high-dimensional 
term, that is by the numerators in Eq. (|22|) . The sign 
of the vertex correction is not affected by the type of 
randomness in the distribution of the scattering impuri- 
ties. Coulomb interaction in the pure Falicov-Kimaball 
model has the same effect on the electrical conductivity 
as variance of fluctuations of the atomic potential in the 
disordered Anderson model, it hinders diffusion. 

We explicitly calculated only the leading high- 
dimensional vertex correction to the Drude zero- 
temperature conductivity in the paramagnetic phase of 



FKM. Quantitatively the correction is almost everywhere 
two orders smaller than the Drude term. It is due to the 
fact that the Drude conductivity diverges in the limit 
of the Fermi gas while the vertex correction asymptot- 
ically approaches a finite value when the interaction is 
switched off. Only close to the metal-insulator transition 
in the electron-hole symmetric case the vertex correction 
is of order of the Drude one, however, they both vanish 
at the transition point. The leading high-dimensional 
contribution to the vertex correction of the mean-field 
conductivity is quantitatively negligible. It is important 
only for determining the sign of the vertex correction. 
To obtain a more realistic values of the vertex correc- 
tions in low-dimensional systems, one has to consider the 
full representation of the vertex correction from Eq. ([^^ 
containing the Cooper pole. Only then we are able to 
include a sizable impact of spatial dimensionality. 

The method for evaluating vertex corrections to the 
mean-field electrical conductivity via an asymptotic ex- 
pansion in high spatial dimensions is universal and is 
suitable for any model with elastically scattered elec- 
trons. We can hence use it to investigate whether the 
Coulomb interaction in FKM can lead to suppression 
of diffusion in low-dimensional systems as well as to 
study an interplay between the Anderson localization 
and Mott-Hubbard metal-insulator transition in the dis- 
ordered Falicov-Kimball model. 
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